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Abstract 

Measurements of primary charged hadron multiplicity distributions are presented for 
non-single-diffractive events in proton-proton collisions at centre-of-mass energies of 
a/s = 0.9, 2.36, and 7TeV, in five pseudorapidity ranges from \rj\ < 0.5 to \tj\ < 2.4. 
The data were collected with the minimum-bias trigger of the CMS experiment dur- 
ing the LHC commissioning runs in 2009 and the 7TeV run in 2010. The multiplicity 
distribution at y^s = 0.9 TeV is in agreement with previous measurements. At higher 
energies the increase of the mean multiplicity with \/s is underestimated by most 
event generators. The average transverse momentum as a function of the multiplic- 
ity is also presented. The measurement of higher-order moments of the multiplicity 
distribution confirms the violation of Koba-Nielsen-Olesen scaling that has been ob- 
served at lower energies. 
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1 Introduction 

The charged hadron multiplicity, or number of primary charged hadrons, n, is a basic global 
observable characterising final states in high-energy-collision processes. The multiplicity dis- 
tribution, P n , is the probability to produce n charged hadrons in an event, either in full phase 
space or in restricted phase space domains. In this paper we report measurements of P„ in non- 
single-diffractive [1J proton-proton collisions, at centre-of-mass energies \/s = 0.9, 2.36, and 
7TeV at the Large Hadron Collider (LHC) [2 J. The measurements are based on events recorded 
by the Compact Muon Solenoid (CMS) [3j experiment, using a minimum-bias trigger. 

Energy-momentum and charge conservation significantly influence the multiplicity distribu- 
tion for the full phase space. The distribution in restricted phase space, which is less affected by 
such constraints, is expected to be a more sensitive probe of the underlying dynamics and can 
be used to better constrain phenomenological models. Comprehensive reviews on the subject 
can be found in HUSH?). The measurements described in this paper are performed for inter- 
vals of increasing extent in pseudorapidity from \ t]\ < 0.5 up to \t]\ < 2.4, where rj is defined as 
— ln[tan(0/2)], and 9 is the polar angle of the particle with respect to the counterclockwise beam 
direction. In these measurements primary charged hadrons are defined as all charged hadrons 
produced in the interaction, including the products of strong and electromagnetic decays, but 
excluding products of weak decays and hadrons originating from secondary interactions. 

Independent emission of single particles yields a Poissonian P n . Deviations from this shape, 
therefore, reveal correlations. These correlations are predominantly short range in rapidity, 
attributed to cluster decays, and reflect local conservation of quantum numbers in the hadroni- 
sation process. In hadron-hadron interactions, additional large long-range rapidity correlations 
are observed, whose magnitude increases with ^Js |6|. In contrast, in e + e~ annihilations such 
long-range correlations are much weaker and practically absent in two- and three-jet event 
samples |7|. 

The mean of the multiplicity distribution, (n), is equal to the integral of the inclusive single- 
particle density in the considered phase-space domain. Higher-order moments of P n measure 
event-to-event multiplicity fluctuations. They are related to the two-particle and higher-order 
inclusive density correlations and provide more detailed dynamical information than that con- 
tained in single-particle inclusive spectra |8l- [T2| . The average transverse momentum of the 
charged particles, (pj), exhibits a positive correlation with the event multiplicity in hadron- 
hadron collisions ([13-19] and references therein). 

Traditionally, the s dependence of P„ and its moments has been much discussed (H [4[ in 
relation to Koba-Nielsen-Olesen (KNO) scaling 11201 121| . In this framework, one studies the 
KNO function Y(z) = (n)P n , where z = n/(n). If KNO scaling holds, Y(z) and the normalised 
moments Cn = (n) 17 are independent of s. 

Throughout this paper, we compare the data with existing measurements at similar or lower 
centre-of-mass energies and with predictions of the multiplicity distribution and its mean value 
from analytical and event generator Monte Carlo (MC) models. The models are based on the as- 
sumption that hadrons are produced via the fragmentation of colour strings. This comparison 
should allow for a better tuning of the existing MC models to accurately simulate minimum- 
bias events and underlying-event effects. 

The next section gives a short description of the CMS detector. Section 3 describes the MC 
models used in the analysis, while Section 4 presents the data samples. The track reconstruction 
and acceptance are explained in Section 5. Section 6 describes the corrections applied to the 
data. Section 7 lists all relevant systematic uncertainties. The results are discussed in Section 8. 
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3 Models 



2 The CMS detector 

A complete description of the CMS detector can be found in [3[. The CMS experiment uses 
a right-handed coordinate system, with the origin at the nominal interaction point (IP), the x 
axis pointing to the centre of the LHC ring, the y axis pointing up, and the z axis along the 
counterclockwise beam direction. 

The central feature of the CMS detector is a superconducting solenoid of 6 m internal diameter, 
providing an axial magnetic field with nominal strength of 3.8 T. Immersed in the magnetic 
field are the pixel tracker, the silicon-strip tracker (SST), the lead tungstate electromagnetic cal- 
orimeter, the brass/ scintillator hadron calorimeter and the muon detection system. In addition 
to the barrel and endcap calorimeters, the steel /quartz-fibre forward calorimeter (HF) covers 
the region 2.9 < \tj\ < 5.2. 

Two of the CMS subdetectors acting as LHC beam monitors, the Beam Scintillation Counters 
(BSC) and the Beam Pick-up Timing for the experiments (BPTX) devices, were used to trig- 
ger the detector readout. The BSCs are located along the beam line on each side of the IP at a 
distance of 10.86 m and are sensitive in the range 3.23 < \tj\ < 4.65. The two BPTX devices, 
which are located inside the beam pipe at distances of 175 m from the IP, are designed to pro- 
vide precise information on the bunch structure and timing of the incoming beams, with a time 
resolution better than 0.2 ns. 

The tracking detector consists of 1440 silicon-pixel and 15148 silicon-strip detector modules. 
The barrel part consists of 3 (10) layers of pixel (SST) modules around the IP at distances rang- 
ing from 4.4 cm to 1.1 m. Five out of the 10 strip layers are double sided and provide additional 
z coordinate measurements. The two endcaps consist of 2 (12) disks of pixel (SST) modules that 
extend the pseudorapidity acceptance to \tj\ < 2.5. The tracker provides an impact parameter 
resolution of about 100 }im and a pj resolution of about 0.7% for 1 GeV/c charged particles at 
normal incidence 11221 1231 . 

3 Models 

The pythia 6 [24 j generator and its fragmentation model tuned to CDF data [25. 26 1, hereafter 
called pythia d6t, is used as a baseline model to simulate inelastic pp collisions. However, 
at 7 TeV a dedicated pythia tune [27] describing better the high multiplicities is used for cor- 
recting the data. Alternative tunings that differ mainly in the modelling of multiple parton 
interactions have also been considered |26l l28l l29| . PHOJET |30l l3"Tfl is used as an alternative 
event generator that differs mainly in the underlying dynamical model for particle production. 
While pythia contains at least one hard scatter per event, particle production in PHOJET, which 
is based on the dual-parton model [32j, is predominantly soft and contains in general multiple- 
string configurations derived from the dual-parton model with multi-Pomeron exchanges [32 J. 
Each Pomeron exchange gives rise to two strings stretched between either valence or sea par- 
tons. At low energies the dominant process is single-Pomeron exchange, which leads to two 
strings stretched between valence quarks and diquarks. With increasing energy, additional 
Pomeron exchanges occur, forming strings stretched between sea partons. Because the sea par- 
tons carry on average only a small fraction of the momentum of the incident hadrons, these 
strings are concentrated in the central rapidity region. These extra strings [33 . 34] are needed to 
explain the KNO scaling violations observed at high energies [35ll36l| , the increase of the cen- 
tral particle density with increasing energy |35ti37| , the (pj) versus n dependence |38ll39| , and 
long-range rapidity correlations [40] . The distribution of the number of exchanged Pomerons 
can be obtained from perturbative Reggeon calculus and unitarity by means of fits to the mea- 
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sured total, elastic, and diffractive cross sections as described in |4lti43| . Other alternatives 
based on the Colour Glass Condensate (CGC) picture with saturated gluons |44| also make 
predictions for the multiplicity dependence of (pj) and for the long-range rapidity correla- 
tions. 

We also compare our measurements with a new fragmentation model implemented in PYTHIA 
8 1 45 1. In addition to having px _ ordered parton showers, rather than an ordering by virtuality it 
differs from its predecessor in the modelling of multiple-parton interactions and the treatment 
of beam remnants and diffraction. The detailed MC simulation of the CMS detector response is 
based on GEANT4 [46] . Simulated events were processed and reconstructed in the same manner 
as collision data. 

4 Data sample 

Near the end of 2009, the CMS experiment collected two datasets of proton-proton collisions at 
centre-of-mass energies of 0.9 and 2.36 TeV. In March 2010 a new running period at y/s = 7 TeV 
started, of which data collected in the first days have been analysed for this paper. The corre- 
sponding inelastic proton-proton interaction rates were about 11, 3, and 50 Hz, respectively, for 
these datasets. At these rates, the fraction of bunch crossings in which two or more minimum- 
bias collisions occurred is negligible |[47l[48l . 

Diffraction is commonly characterised by one (single-diffraction) or two (double-diffraction) 
colourless exchanges resulting in the observation of a large rapidity interval devoid of any 
hadron activity (rapidity gap). All results presented in this paper refer to inelastic non-single- 
diffractive (NSD) interactions and are based on an event selection that retains a large frac- 
tion of the non-diffractive (ND) and double-diffractive (DD) events, while disfavouring single- 
diffractive (SD) events. 

The trigger and offline event selection are nearly identical to those used in |47ll48l ; about three 
times more events were used in this analysis at y/s = 0.9 and 7 TeV. The trigger required a 
signal in any of the BSC scintillator counters, in coincidence with either of the two BPTX devices 
indicating the presence of at least one proton bunch crossing the IP. Beam halo backgrounds are 
reduced by using the timing information from the BSC counters at opposite ends of the CMS 
detector. Additional beam-induced backgrounds are removed by requiring the cluster sizes in 
the pixel detector to be consistent with a single primary vertex, as described in [47] . NSD events 
are selected by requiring at least one HF calorimeter tower with more than 3 GeV of total energy 
in each of the positive-z and negative-z HF calorimeters. Finally, a reconstructed primary vertex 
is required with the z coordinate within ±15 cm of the centre of the beam collision region. 

For the final analysis, totals of about 132, 12, and 442 thousand events were retained at y/s = 
0.9, 2.36, and 7 TeV, respectively. Event yields at various stages of the selection are shown in 
Tabled] 

5 Track reconstruction and acceptance definition 

The barrel and endcap pixel and SST detectors are used in the reconstruction of tracks within 
an acceptance of \t]\ < 2.5. Due to a large drop in reconstruction efficiency near the limits of 
this range, we restrict the computation of the multiplicity spectra to the region | rj | < 2.4. 

In proton-proton collisions at the LHC, the events selected by minimum-bias triggers involve 
predominantly soft interactions and contain mostly particles with small transverse momenta. 
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5 Track reconstruction and acceptance definition 



Table 1: Event yields in each data sample after sequential trigger and event selection. 



Selection 


y/s (TeV) 


0.9 


2.36 


7 


beam background rejection + LI trigger 
all preceding + forward calorimeters 
all preceding + primary vertex 


254666 
146658 
132294 


18739 
12019 
11674 


610549 
500077 
441924 



These are reconstructed by extending the standard tracking algorithms of the CMS experiment, 
which are based on a combinatorial track finder Ij22l that performs multiple iterations. Hits that 
can be assigned unambiguously to tracks in one iteration are removed from the collection of 
tracker hits to create a smaller collection that can be used in the subsequent iteration. This iter- 
ative procedure was further optimised for primary track reconstruction with pj > 100 MeV/c 
in minimum-bias events [47]. The reconstruction efficiency estimated by means of the de- 
tector simulation, exceeds 90% for tracks with pj > 500MeV/c and drops below 70% for 
tracks with pj < lOOMeV/c. Contamination from mis-reconstructed tracks is below 5% for 
pj < 500MeV/c. 

After three iterations of the combinatorial track finder, the position of the primary vertex is 
recomputed and then used as an additional constraint in a refit of all previously reconstructed 
tracks, thus improving the overall resolution in t] and pj. An agglomerative clustering al- 
gorithm followed by a Gaussian mixture model [|49| were applied in order to optimise the 
vertex-finding efficiency. 

The contamination due to decays of long-lived particles (denoted as V° decays) is dominated 
by charged pions and protons originating from Kg and A decays. The Kg production rate is 
roughly 5% of that for all charged particles |50l[ , and the resulting charged secondaries amount 
to 7% of all 7T ± . The A production is measured to be 43% of Kg [51, 52 j in this kinematic do- 
main, yielding another 3% of charged secondaries consisting of protons and pions. In order 
to reduce the contamination of these V° decay products and secondaries produced in interac- 
tions of charged particles with the material of the detector, we require all reconstructed tracks 
to be associated with the primary vertex. This is done by selecting tracks with a small impact 
parameter with respect to the primary vertex position both in the transverse plane and along 
the z axis [47J. Only 1.5% of Kg and A decay products pass these selections, resulting in a final 
contamination of 0.2%. 

The K ± / 7r ± ratio is known to be fairly constant over a wide range of centre-of-mass energies 
and is between 8 and 12% [|50 | . The K ± have a lower reconstruction efficiency at low pj than 
pions and mismodeling of the K ± / 7T ± ratio could result in a change in the multiplicity dis- 
tribution. However a doubling of the K ± / zr ± ratio yields a negligible shift of 0.25% in the 
multiplicity average. This is expected because the (pj) of K ± is substantially higher than n , 
and K ± therefore contribute very little in the pj range where the reconstruction efficiencies dif- 
fer. Finally, only tracks with a relative uncertainty on their measured transverse momentum 
smaller than 10% are selected; this requirement rejects mainly low-quality and badly recon- 
structed tracks. 
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Table 2: The percentage of single-diffractive events (SD/ inelastic) at each centre-of-mass energy 
as predicted by pythia d6t and PHOJET. 



Generator 


y/s (TeV) 


0.9 


2.36 


7 


PYTHIA D6T 
PHOJET 


22.5 
18.9 


21.0 
16.2 


19.2 
13.8 



6 Corrections 

Due to the requirement of significant activity in both ends of the HF, the event selection ac- 
ceptance for SD events is small: 5% at 0.9 TeV and 7% at 7 TeV. This contribution is therefore 
subtracted based on simulated pythia events. The pythia and PHOJET predictions of the ini- 
tial fractions of SD events differ substantially, as seen from Table [2] The difference of the two 
predictions is taken as the systematic uncertainty related to the SD subtraction. 

It is customary to normalise the charged hadron multiplicity distribution P n = a n / a, where 
a n denotes the cross section for a fixed multiplicity n, to either the total inelastic cross section 
or the NSD cross section. For the results presented in this paper the normalisation factor a 
corresponds to the latter. 

The minimum-bias trigger and NSD selection unavoidably introduce a bias in the measured 
charged hadron multiplicity. Furthermore, a fraction of the events are removed by the require- 
ment of a good quality primary vertex. These effects result in an accepted multiplicity given 

Tn = €n ' Pn/ (1) 

where e„ is the trigger and event reconstruction efficiency for multiplicity n. This efficiency is 
close to 100% for multiplicities larger than n = 20 and drops gradually to 40% for n = 1. 

Due to inefficiencies in track reconstruction and acceptance, the creation of secondary particles 
by the interaction of primaries with the beam pipe and the detector material, and the presence 
of decay products of long-lived hadrons, one will in general not measure the true multiplicity 
n but a statistically related quantity m. The statistical distribution O m of this observed multi- 
plicity is related to the true accepted multiplicity distribution by the linear relation 

Om = ^^Rm,n ' T n . (2) 

n 

The problem of inverting the response matrix R m ,n/ here taken from MC, is well known and 
extensive literature on the topic exists 1153 [ [54| . When dealing with limited event samples, 
an algebraic inversion of R m/tl turns out to be impossible and leads to unstable results. We 
therefore use a Bayesian unfolding method, as described in [55|. The unfolding procedure in- 
troduces large statistical correlations between adjacent bins of the multiplicity spectrum. The 
full covariance matrix of the unfolded multiplicity spectrum was calculated using a resampling 
technique fl54| . 

The average transverse momentum (px) as a function of n was measured as well. For each bin 
in raw multiplicity, a Monte Carlo based correction factor (p| en )/ (Px c ) was applied to convert 
the measured (px) to the corresponding value for primary charged hadrons. Subsequently, the 
response matrix R m/n , weighted by the corrected data, was applied to correct the multiplicities. 

The correctness of the unfolding procedure was verified using MC events generated with our 
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7 Systematic uncertainties 



baseline pythia d6t tune where we compared the multiplicity distribution of generated pri- 
mary hadrons with that of reconstructed tracks after unfolding. At all energies considered, 
an excellent agreement with the generated hadrons is achieved, proving the stability of the 
procedure. 

The track reconstruction efficiency of the minimum-bias tracking drops drastically for pj < 
lOOMeV/c, while the mis-reconstructed track rate increases. Rather than correcting for these 
effects using MC simulations, we reconstruct the pj spectrum in data themselves and calculate 
the fraction of charged hadrons with pj < 100 MeV/c by extrapolating the measured spectrum 
in data, using a parametrisation based on an exponential of a third-degree polynomial in pj. 
The fraction of charged hadrons added by this correction ranges between 5% and 7% depend- 
ing on the centre-of-mass energy and the pseudorapidity interval under study. The functional 
form used to extrapolate to the lowest transverse momenta introduces an uncertainty of 1% 
on this fraction and is taken into account in the systematic uncertainties. The effect of all the 
correction procedures on the measured raw multiplicity distribution is illustrated in Fig. [T] for 
a centre-of-mass energy of 7 TeV. 




Figure 1: A comparison of the uncorrected and fully corrected multiplicity distribution at 
s/s = 7 TeV for \t]\ < 2.4 . The uncertainties before corrections are statistical only, while af- 
ter corrections the statistical and systematic uncertainties are added in quadrature. 



7 Systematic uncertainties 

Four main sources of systematic uncertainty contribute to the total uncertainty on the multi- 
plicity distribution P n and (pj) versus n: the uncertainty on the trigger and event selection 
efficiency, the uncertainty on the tracking efficiency and acceptance, the model dependence of 
the response matrix, and the model dependent SD subtraction. All four are discussed below. 
The effects of misalignment of the tracking detector, dead sensors, and the uncertainty on the 
vertex position are much smaller and are contained within the overall tracking systematic un- 
certainty. For (pj) versus n, the effect of the event selection and SD subtraction was observed 
to be negligible. 
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Table 3: Summary of systematic uncertainties on the track reconstruction. 



Source 


tracking uncertainty (%) 


Tracking efficiency 


2.0 


Acceptance 


1.0 


Pixel hit efficiency 


0.3 


Pixel cluster splitting 


0.2 


Correction for secondaries 


1.0 


Misalignment 


0.1 


Beam halo 


0.1 


Multiple track counting 


0.1 


Mis-reconstructed tracks 


0.5 


Total 


2.5 



Trigger and event selection efficiency. The corrections for the trigger and event selection effi- 
ciency are based entirely on MC simulation. The largest impact on the overall efficiency is that 
of the HF coincidence requirement. A cross-check of the multiplicity dependent efficiency with 
zero-bias events, which are by definition not biased at all (random trigger on collisions), shows 
good agreement within statistical errors between data and MC. A relative shift of the efficiency 
correction factors by for low multiplicities, decreasing to ±1% for n ~ 20, covers the trigger 
efficiency measured in zero-bias data. This leads to a maximum 5% systematic uncertainty. 

Tracking efficiency and acceptance. A correct description of the tracking efficiency in the MC 
simulation of the detector is essential for obtaining a correct response matrix. At low trans- 
verse momenta, the efficiency drops due to a loss of hits on the tracks that are stopped within 
the tracking volume. As in ||47| , we assign a 2% uncertainty on this efficiency. The remaining 
contamination is mostly due to secondary tracks originating from interactions with the mate- 
rial of the LHC beam pipe around the interaction point. This is estimated in [47j to be no more 
than 1%, and is confirmed based on [56 j. These uncertainties, together with smaller contribu- 
tions from misalignment of the tracking detector, beam-halo background, multiple counting of 
tracks, and mis-reconstructed tracks are added in quadrature to produce a total tracking uncer- 
tainty of 2.5%. All correction factors related to track reconstruction are summarised in Table [3] 
As was discussed in Section [5J the contamination from V° decays after associating tracks with 
the primary vertex is small (0.2%) and already included in the systematic uncertainty for sec- 
ondary tracks. The difference in reconstruction efficiency for charged kaons and pions also has 
a negligible impact on the measured multiplicity distribution. Finally, the extrapolation uncer- 
tainty from px = lOOMeV/c to zero is 1% (Section [6j. For (pj) versus n, this was taken into 
account by changing the mean transverse momentum by the average pj of tracks lost during 
the event reconstruction folded with the tracking efficiency uncertainty. The overall effect of 
this systematic uncertainty is less than 10% for small multiplicities, but can reach up to 30% for 
the high multiplicities, where it is the main source of systematic uncertainty. 

Model dependence. The baseline MC model that is used to unfold the multiplicity distribution 
underestimates the single-particle densities at zero rapidity by 10% relatively to \t}\ ~ 1.5 at 
yjs = 0.9 and 2.36 TeV, but correctly describes the 7TeV data shapes. This discrepancy affects 
the response matrix used in the unfolding procedure. The effect on the multiplicity distribution 
is estimated to be at most 3% by means of an alternative MC tune. The robustness of the 
unfolding procedure was verified by unfolding pseudo-data generated with PHOJET using a 
response matrix constructed with PYTHIA. The induced variation on P n is below 3% and is 
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8 Results 



therefore considered to be covered already by the systematic uncertainty related to modelling 
of the single-particle densities. 

SD subtraction. Finally the modelling of SD events was found to be different in pythia and 
PHOJET, both in the predicted relative fraction and in the predicted multiplicity distribution. 
The largest variation is obtained by subtracting the PHOJET multiplicity distribution for SD 
events, using the event fraction predicted by pythia and vice versa. The impact can be as 
large as 20% in some low-multiplicity bins. At high multiplicities this uncertainty is below 5%. 

All systematic uncertainties are calculated in each multiplicity bin for upward and downward 
effects on the multiplicity distribution and added in quadrature to yield a total systematic un- 
certainty. The total uncertainty on P n is below 10% for a large part of the multiplicity distribu- 
tion, but increases to 40% at the lowest and highest multiplicities. 

8 Results 

8.1 Charged hadron multiplicity distributions 

The NSD charged hadron multiplicity distributions are measured in increasing ranges of pseu- 
dorapidity from \tj\ < 0.5 to \r]\ < 2.4. The fully corrected results at \/s = 0.9, 2.36, and 
7 TeV are compared in Fig. [2] with earlier measurements in the same pseudorapidity ranges 
performed by the UA5 I|35l l36 | and ALICE Ii57l l58f collaborations. Our measurements were 
also compared with results obtained with a CMS cross-check analysis (not shown) on data at 
a/s = 0.9 and 7 TeV, using a tracklet-based tracking algorithm as in Ref. [l4"7| . With a recon- 
struction efficiency exceeding 90% for pj > 50 MeV/ c, the latter provided a cross-check of the 
extrapolation for tracks below pj < lOOMeV/c, including the use of data without magnetic 
field at 7 TeV. All measurements agree well within their total uncertainties. 

In the largest pseudorapidity interval of 1 77 1 < 2.4, there is a change of slope in P n for n > 20, 
indicating a multicomponent structure as was discussed in [59| and [60 1 in terms of multiple- 
soft-Pomeron exchanges. This feature becomes more pronounced with increasing centre-of- 
mass energies, notably at a/s = 7 TeV. 

An extensive range of tunes [l26l - l2"9"| based on the pythia 6 fragmentation model have been de- 
veloped. They differ mainly in their parametrisation of the multiple-parton-interaction model. 
Some reproduce the charged hadron multiplicities better than others, but none is able to give 
a good description simultaneously at all the centre-of-mass energies and in all pseudorapidity 
ranges. For clarity, only the baseline tune D6T is shown in comparison with other models hav- 
ing a different physical description of soft-particle production such as PHOJET 11301 |3T|| and the 
new fragmentation model of PYTHIA 8 |45| . 

A comparison of our measurements with these three classes of models is shown in Fig. [3] for 
all charged hadrons and for those with pj > 500MeV/c. pythia d6t underestimates drasti- 
cally the multiplicity at all measured energies but improves when pj > 500MeV/c is required. 
pythia 8 is the only model that gives a reasonable description of the multiplicity distribution 
at all energies, but tends to overestimate the multiplicity at 7 TeV when pj > 500 MeV/ c is 
required. PHOJET produces too few charged hadrons overall but gives a good description of 
the average transverse momentum (pj) at fixed multiplicity n, as illustrated in Fig. |4] Among 
the three classes of models, pythia 8 gives the best overall description of the multiplicity dis- 
tribution and the dependence of the average transverse momentum on n. Inspired by |6l| we 
fit a first-degree polynomial in ^Jn to the multiplicity dependence of (pj) (n) for n > 15 at each 
energy, yielding a good description which is valid at all three energies (Fig. B. The ratios of 
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Figure 2: The fully corrected charged hadron multiplicity spectrum for \rj\ < 0.5, 1.0, 1.5, 2.0, 
and 2.4, (a) at y/s = 0.9 TeV, (b) 2.36 TeV, and (c) 7TeV, compared with other measurements in 
the same rj interval and at the same centre-of-mass energy |35ll36ll57ll58| . For clarity, results in 
different pseudorapidity intervals are scaled by powers of 10 as given in the plots. The error 
bars are the statistical and systematic uncertainties added in quadrature. 



the data obtained at 7 and 2.36 TeV with respect to the data at 0.9 TeV show that the rise of the 
average transverse momentum with the multiplicity is roughly energy-independent . 

All previous observations seem to indicate that the Monte Carlo models produce too few parti- 
cles with low transverse momenta, especially at 7 TeV. The pythia models tend to compensate 
for this by producing too many particles with high transverse momentum, which is related to 
the modelling of semi-hard multiple-parton interactions. 



10 



8 Results 




Figure 3: The charged hadron multiplicity distributions with \t]\ < 2.4 for (a) pj > and (b) 
p T > 500 MeV/c at y/s = 0.9, 2.36, and 7 TeV, compared to two different pythia models and the 
PHOJET model. For clarity, results for different centre-of-mass energies are scaled by powers of 
10 as given in the plots. 
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Figure 4: (a) A comparison of (pj) versus n for < 2.4 with two different PYTHIA models 
and the PHOJET model at a/ s = 0.9, 2.36, and 7 TeV. For clarity, results for different energies 
are shifted by the values of a shown in the plots. Fits to the high-multiplicity part (n > 15) 
with a linear form in ^Jn are superimposed, (b) The ratios of the higher-energy data to the fit 
at 0.9 TeV indicate the approximate energy independence of (pj) at fixed n. 



8.2 Violation of KNO scaling 
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8.2 Violation of KNO scaling 

The multiplicity distributions are shown in KNO form in Fig. [5] for a large pseudorapidity 
interval of\t]\ < 2.4, where we observe a strong violation of KNO scaling between = 0.9 TeV 
and 7 TeV, and for a small pseudorapidity interval of \t]\ < 0.5, where KNO scaling holds. 
Scaling is a characteristic property of the multiplicity distribution in cascade processes of a 
single jet with self-similar branchings and fixed coupling constant II62H69II . 




z = n/<n> z = n/<n> 

Figure 5: The charged hadron multiplicity distributions in KNO form at \Js = 0.9 and 7 TeV in 
two pseudorapidity intervals, (a) \t]\ < 2.4 and (b) \tj\ < 0.5. 

The validity of KNO scaling is shown more quantitatively in Fig. |6]by the normalised order-^ 
moments Cq of the multiplicity distribution, complemented with measurements at lower ener- 
gies |70Tl72| . For \r/\ < 2.4 the values of C q increase linearly with logs, while for \t]\ < 0.5 they 
remain constant up to q = 4 over the full centre-of-mass energy range, as illustrated by the fits 
in Fig. [6] 

Multiplicity distributions for e + e annihilations up to the highest LEP energies show clear 
evidence for multiplicity scaling, both in small ranges (At] < 0.5), in single hemispheres, and 
in full phase-space. However, at LEP energies, scaling is broken for intermediate-size ranges 
where, besides two-jet events, multi-jet events contribute most prominently 1 73] - t77| . 

For hadron-hadron collisions, approximate KNO scaling holds up to ISR energies l78l|79|, but 
clear scaling violations become manifest above a/s ~ 200 GeV both for the multiplicity distri- 
butions in full phase space and in central pseudorapidity ranges Il59l [70l l80ll8TH . In pp colli- 
sions, and for large rapidity ranges, the UA5 experiment was the first to observe a larger than 
expected high-multiplicity tail and a change of slope ||59l l72| , which was interpreted as evi- 
dence for a multi-component structure of the final states f34ll60l[8 2]|. Our observation of strong 
KNO scaling violations at \/s = 7 TeV, as well as a change of slope in P n , confirm these earlier 
measurements. 

All these observations, together with the sizable growth with energy of the non-diffractive 
inelastic cross section, point to the increasing importance of multiple hard, semi-hard, and soft 
partonic subprocesses in high energy hadron-hadron inelastic collisions I6l l32ll34ll59ll83ll84ll . 
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Figure 6: Fits of the log s dependence of the normalised moments Cn of the multiplicity distri- 
bution for (a) \y\ I < 2.4 (assuming linear dependence) and (b) \r\\ < 0.5 (assuming no depen- 
dence), including data from lower energy experiments II70H72|| . For \/s = 0.9 TeV, data from 
experiments other than CMS were drawn shifted to lower a/s for clarity. 



8.3 Energy dependence of the mean multiplicity 

The mean multiplicity (n) is the first moment of the multiplicity distribution and is equal to 
the integral of the corresponding single-particle inclusive density in the rj interval considered. 
The mean multiplicity (n) is observed to rise with increasing centre-of-mass energy in hadron- 
hadron collisions I35l- l37ll70l[78ll79ll85ll86| . The same behaviour is also observed in e + e colli- 
sions, in deep-inelastic scattering [87J, and in heavy ion collisions IQX 

Our measured mean multiplicity is compared with experimental data obtained at lower ener- 
gies in Fig. [7] Recent Regge-inspired models |4iTl43t predict a power-like behaviour among 
which only Ref. [42J describes the highest energy data very well. Parton saturation models 
(such as [44j) predict a strong rise of the central rapidity plateau as well. Table [4] gives an 
overview of (n) for the data and for the pythia d6t, pythia 8, and PHOJET event generators. 



As in Section 8.1 pythia d6t produces on average too few particles per event at all three en- 
ergies. PHOJET is consistant with the data within uncertainties for \/s = 0.9 TeV, but is not able 
to predict properly the mean multiplicity at higher energies, pythia 8 describes best the 7 TeV 
data, but underestimates (n) systematically at all energies. 
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Table 4: Mean multiplicity for data, pythia d6t, pythia 8, and PHOJET for \rj\ < 2.4 at each 
centre-of-mass energy. For data, the quoted uncertainties are first statistical, then upward and 
downward systematic. 
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Figure 7: The evolution of the mean charge multiplicity with the centre-of-mass energy for 
\t] | < 2.4, including data from lower-energy experiments for \t]\ < 2.5 l37l |70T[72| . The data 
are compared with predictions from three analytical Regge-inspired models [41 -43 J and from 
a saturation model 



9 Conclusions 

The charged hadron multiplicity distributions of non-single-diffractive events were measured 
from an analysis of the minimum-bias datasets collected by CMS at three centre-of-mass ener- 
gies: \/s = 0.9, 2.36, and 7 TeV. The excellent tracking capabilities of the silicon pixel and strip 
detectors of CMS, combined with an optimised tracking and vertexing algorithm, allow the 
reconstruction of charged tracks down to = lOOMeV/c with high efficiency and low back- 
ground contamination. A full correction for detector resolution and acceptance effects and an 
extrapolation to zero transverse momentum yield measurements of the charged hadron mul- 
tiplicity distribution for increasing central pseudorapidity ranges from \rj\ < 0.5 to \t]\ < 2.4, 
which can be compared with models of soft-particle production and with experimental data at 
lower energies. 

Although some event generators provide an adequate description of Tevatron and LEP data, 



14 



9 Conclusions 



none is able to describe simultaneously the multiplicity distributions and the px spectrum at 
y/s = 7 TeV. In general, models predict too few low-momentum particles, indicating that by 
increasing the amount of multiple-parton interactions one effectively introduces too many hard 
scatters in the event. 

The change of slope in P n in the widest central pseudorapidity intervals observed at \/s = 
7 TeV, combined with the strong linear increase of the C q moments, indicates a clear violation 
of KNO scaling with respect to lower energies. This observation merits further studies. 
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